Level statistics of a pseudo-Hermitian Dicke model 
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A non-Hermitian operator that is related to its adjoint through a similarity transformation is 
defined as a pseudo-Hermitian operator. We study the level-statistics of a pseudo-Hermitian Dicke 
Hamiltonian that undergoes Quantum Phase Transition (QPT). We find that the level-spacing 
distribution of this Hamiltonian near the integrable limit is close to Poisson distribution, while it is 
Wigner distribution for the ranges of the parameters for which the Hamiltonian is non-integrable. 
We show that the assertion in the context of the standard Dicke model that QPT is a precursor to 
a change in the level statistics is not valid in general. 
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The study on statistics of energy-levels in quantum 
many-body systems has a long history [l|, Q . The statis- 
tical analysis based on Random Matrix Theory (RMT) 
has been apphed to characterize quantum chaos and 
to investigate the integrabihty of a quantum system. 
In particular, it has been conjectured [3] that the level- 
spacing distribution of an integrable Hermitian Hamil- 
tonian should be described by the Poisson distribution: 
Pp(s) — exp(— s). On the other hand, if the system is 
non-integrable, the level-spacing distribution of the Her- 
mitian Hamiltonian should be given by the Wigner distri- 
bution, i.e., the Wigner surmise for the Gaussian Orthog- 
onal Ensemble (GOE): Pw(s) = (7rs/2) exp(-7rsV4)- 
Although there is no rigorous proof of the Bohigas- 
Giannoni-Schimdt(BGS) conjecture^ for quantum sys- 
tems, it has been numerically confirmed for a variety of 
many-body Hamiltonian and also theoretically in the 
semi-classical limit [3]. 

The RMT without the constraint of Hermiticity was 
introduced by Ginibre Q and currently, is an active 
field of research fo*]. The non-Hermitian RMT exhibits 
generic statistical behavior of quantized dissipative sys- 
tems. The integrable case corresponds to the Poisson 
process on the plane, while a cubic repulsion is a sig- 
nature of quantum chaotic scattering [6[. An interest- 
ing result due to Ginibre Q is that the probability den- 
sity function for the eigenvalues of real Gaussian random 
non-symmetric matrices with all the eigenvalues being 
real is identical to the GOE and consequently, the level- 
spacing distribution is given by P^{s). If all the eigen- 
values of a real non-symmetric matrix M are real, it can 
be shown that the same matrix can be mapped to its 
transpose through a similarity transformation. In par- 
ticular, = {XX'^)-'^M{XX'^), where the real ma- 



* Electronic address: 'deguchi@phys.ocha.ac.jpl 
^Electronic address: pijushkanti.ghosh@visva- bharati.ac.in| 
■f Electronic address: |kudo.kazue@ocha.ac.jp| 



trix X diagonalizes M with entirely real eigenvalues E, 
i.e., M = XEX^^. This shows that Ginibre's ensem- 
ble of real non-symmetric matrices belong to the class of 
operators known as pseudo-Hermitian operator, i.e., an 
operator that is related to its adjoint through a similarity 
transformation. 

The study on pseudo-Hermitian operators have re- 
ceived considerable attention recently in connection with 
the pioneering work of Bender and Boettcher f?\ show- 
ing that non-Hermitian operators with unbroken VT- 
symmetry admit entirely real spectra. This has opened 
up several new directions [1, [HI, in the study of 
pseudo-Hermitian operators. One of the significant de- 
velopments is the construction of pseudo-Hermitian RMT 
with pseudo- unitary symmetry [l3| . The degree of level 
repulsion is different from that of previously known Gaus- 
sian orthogonal, unitary and symplectic ensembles and it 
seems to point out a new universality class. Moreover, 
unlike the Ginibre ensembles (except for the exceptional 
case discussed above) , non-Hermitian RMT with pseudo- 
unitary symmetry describes non-dissipative systems. 

In spite of all the above developments, a criterion 
to characterize quantum chaos and to investigate the 
integrabihty of a pseudo-Hermitian operator using the 
level-statistics based on RMT is still lacking. Pseudo- 
Hermitian operators with entirely real spectra can be 
shown to be Hermitian with respect to some modified 
inner product in the Hilbert space Q. The effect of 
the modified inner product in the Hilbert space is to 
have a modified symplectic structure for the correspond- 
ing classical system. It may be noted here that a fixed 
modified inner product in the Hilbert space or the corre- 
sponding symplectic structure is not universal for pseudo- 
Hermitian systems. It varies from one system to another 
and in general, it is a difficult problem to identify the 
proper inner product for a given pseudo-Hermitian sys- 
tem admitting an entirely real spectrum. It is thus not 
clear a priory whether the standard semi-classical analy- 
sis in support of the BGS conjecture for quantum systems 
with standard inner product in the Hilbert space will re- 
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main unchanged for generic pseudo-Hermitian systems 
or not. In the absence of any theoretical support for the 
vahdity of the BGS conjecture for pseudo-Hermitian sys- 
tems, it may be worth fooking for numerical evidences. 

The purpose of this paper is to present numerical evi- 
dence to show that the level-spacing distribution of a non- 
integrable pseudo-Hermitian Dicke Hamiltonian with an 
entirely real spectrum is described by the Wigner distri- 
bution. On the other hand, it approaches to the Poisson 
distribution for the parameters of the model close to the 
integrable limit. 

We consider a pseudo-Hermitian Dicke Hamiltonian 
(DH) that has been shown recently to undergo QPT ^] , 
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where ujq, a, (3, 7, 5, ^1, ^2 are real parameters and j is 
the total spin-angular momentum. The operators a, 
are the standard bosonic annihilation-creation operators 
and Jz, J± are the generators of the SU{2) algebra, 

[a, a)] 1, 

[J+,J_]=.2J„ [Jz,J±] =±J±. (2) 
The Hamiltonian H commutes with the parity operator 

n, 



H 



i-rrN 



N ^ a^a + Jz + j. 



(3) 



The eigenstates of H have definite parity depending on 
whether the eigenvalues of the operator N are odd or 
even. The Hamiltonian ([T]) reduces to the standard DH 
for £,1 = £,2 = and a — f3 = 'y = 5. The DH has been 
studied extensively from the viewpoint of QPT Il4l . [ip . 
[iGj . level-statisti cs fl6| . quantum entanglement jlTl. Il8||. 
exact solvability [l9| and two-dimensional semiconductor 
physics [20j . 

The non-Hermitian Hamiltonian H can be mapped 
to a Hermitian Hamiltonian Ti. = pHp^^ through a 
similarity transformation when the following relation is 
satisfied JJi], 



a (5-/37 = 0. 



(4) 



The operator p and Ti, have the following forms jl2l|: 
p — exp 



a 7 




J+a) 



(5) 



The Hamiltonian H that is non-Hermitian under the 
Dirac-Hermiticity condition becomes Hermitian with re- 
spect to the modified inner-product defined in the Hilbert 



space as, {{u, v))^ 
p^ . In particular. 



[u^rj+v), where the metric 774 



{u\Hv) ^ {Hu\v), {{u\Hv))^^ - {{Hu\v))r 



(6) 



Thus, with the modified inner-product, the results of a 
Hermitian Hamiltonian follow automatically. Note that 
{u\nv) = {nu\v). We refer to Ref. [12] for further de- 
tails. 

In this system, when j is finite, the parity H is a good 
quantum number. Two states with different parity do 
not interact with each other. In other words, we can 
concentrate on the states with either positive or negative 
H. Here, we consider the positive-parity states. 

The level-spacing distributions are given by the prob- 
ability function P{s) of nearest-neighbor spacings Si — 
Xi+i — Xi, where unfolded eigenvalues. In order 

to characterize the level-spacing distribution, we employ 
the quantity. 



Pw{s)]ds 
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where so = 0.4729 ... is the intersection point of Pp{s) 
and Pvv(s)- We have 77 = 1 when P{s) = -Pp(s), and 
ri = when P{s) = Pw(s). 

In the following, we set uj — luq — I and Ci — £2 — for 
convenience. A further choice of 7 = ^ essentially fixes 
6 as S = ^ due to the pseudo-Hermiticity condition ([5]), 
where n{j^ —1) is a real number. With this parametriza- 
tion, the Hamiltonian H can be re-written as, 

H = a'^a+Jz-i — p= ( aJ_a^ -I- pj+a ^ J_a -t- —Jj^a) 

y'^j \ n n 

(8) 

The equivalent Hermitian Hamiltonian Ti has the follow- 
ing form: 



Ti. = a 




n 



(9) 

The total spin-angular momentum j should be large 
enough to obtain proper results of level statistics. If j 
is very small (j ~ 1), because of a kind of finite size ef- 
fects, level statistics shows no universal ensembles [l6j . 
In our numerical calculation, j = 10 unless specifically 
mentioned. 

Figure [1] exhibits the phase diagram of 77 for H in Eq. 
([5]) with n = 1. At the critical line af3 = 1/4, 77 rapidly 
changes. For af] < 1/4, level statistics is almost Poisso- 
nian as seen in Fig. [2ta) [2l[. As a/? increases, P{s) 
changes from the Poisson to Wigner distributions: It 
gives an intermediate distribution, e.g. Fig. ^h). For 
a/3 > 1/4, level-spacing distribution is almost given by 
the Wigner distribution, as shown in Fig. [2l^c). However, 
as a/3 increases further, rj gradually increases. In other 
words, level statistics gradually changes from the Wigner 
distribution to Poissonian one again as a/3 becomes very 
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FIG. 1: (color online) Phase diagram of ri, which is defined 
by Eq. ([7]), in a special case where a = 7 and /3 = S. The 
solid curve corresponds to the critical line a/3 = 1/4. Solid 
circles correspond to the level-spacing distributions in Fig. [21 
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FIG. 2; Level-spacing distributions for the Hamiltonian ^ 
with n = 1. (a) Almost Poisson distribution at a = 7 = 0.1 
and (3 = 5 — 0.2, (b) intermediate distribution at a = 7 = 0.5 
and (3 = 5 = 0.3, (c) Wigner distribution at a = 7 = 0.5 
and (3 = 5 = 0.7, and (d) deformed Wigner distribution at 
a = 7 = 4 and (3 = 5 = 4.5. 



large. Figure [2l^d) is an example of an intermediate dis- 
tribution for large a(3. The behavior of 77 along the line 
a = (3 in Fig. [T] corresponds to that of Ref. [TB] ■ 

The Hamiltonian Ti, in Eq. ([9]) corresponds to the 
standard DH for n = 1 and has been studied in some 
detail in Ref. [l6|. Based on the results of [IB] and the 
numerical findings for H as described above, we suggest 
the following: 

(i) The criteria to distinguish between integrable and 
non-integrable phases of a Hermitian Hamiltonian are 
valid also for a pseudo-Hermitian Hamiltonian. The 
quantum Hamiltonian Ti. and H have the identical 
eigenvalues, since they are related to each other through 
a similarity transformation. Thus, both H and Ti. show 
similar changes in the level-spacing distributions as a 
function of ap. 

(ii) The onset of quantum chaos in a pseudo- Hermitian 
Hamiltonian is manifested by a change in the level 
statistics from Poissonian to Wigner distribution. It 
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FIG. 3: (color online) Phase diagram of rj, which is defined by 
Eq. (0 for 7 = a/n and 5 = (3/n; (a) n = 1/2, (b) n = 2, (c) 
n = 1/3, and (d) n = 3. The solid curve in each panel is given 
by q/3 = n^/(n -I- 1)^. (e) Dependence of rj plotted against 
n for a = 1.25 and (3 = 0.8. The parameter combination is 
shown in (a) - (d) by the points (x). 



is known that the semi-classical Hermitian Hamilto- 
nian corresponding to Ti shows chaotic behavior for 
a(3 > J and regular periodic orbits are obtained for 
aP < J [1^. It is also the case for the Hamiltonian 
H in the semi-classical limit. In fact, a non-Hermitian 
Hamiltonian and its equivalent Hermitian Hamiltonian 
describe the same physics in the classical limit within 
the formalism of pseudo-Hermitian quantum physics and 
the correspondence principle Q. 



We now present numerical results for other values of n. 
Figures ^a.) and ^h) are phase-diagrams of rj for n = ^ 

and n — 2, respectively. Changes in level-statistics ap- 

2 

pear around the critical-line determined by a(3 = ^^^^ 
for both the cases, i.e., the physical picture is identical to 
the case of n = 1. However, the change of rj around the 
critical line is small in Figs.[3]Jc) and[2Id), which exhibit 
phase-diagrams of 77 for n — ^ and n — 3, respectively. In 
fact, level statistics does not show clear Wigner behavior 
for aB > I " with n = ^ and n = 3. 

In order to see the change of 77 clearly, we depict 77 as a 
function of ti for a(3 = 1 in Fig.[3l^e), where the curve of 
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r] for J = 10 is compared with that for j = 30. The non- 
zero r] behavior for n > 2 and n < ^ in Fig.[3]Je) indicates 
that Poisson behavior continues beyond the critical Hne 
in the phase diagram of rj for those regimes of n. The 
behavior is a manifestation of the fact that the system 
becomes close to the integrable limit, if |n| <C 1 or \n\ ^ 
1. One plausible explanation to understand the origin 
of the precise critical values of n (i.e., ^ and 2) may 
lie in the non-applicability of perturbation techniques by 
treating either the counter-rotating ( for 1 < n < 2 ) or 
the rotating terms ( for 1 > rt > i) as perturbation. 

A comment is in order at this point. The pscudo- 
Hermitian Dicke model is known to undergo QPT with 
the critical-line determined by the equation a(3 — 
(^^^^lyi [13] • We suggest that the value of 77, which char- 
acterizes the level statistics, is a possible measure to es- 
timate the onset of the QPT not only for the Hermitian 
Hamiltonian 7i, but, also for the quasi-Hermitian Hamil- 
tonian H. However, the usual assertion (16] within the 
context of the standard Dicke model that QPT is a pre- 
cursor to a change in the level statistics is not valid in 
general for the pseudo-Hermitian Dicke model which has 
a larger parameter space. According to the relevant as- 
sertion in Ref . [3] , level statistics is given by the Wigner 

2 

distribution for af3 > (^n+iy ■ For any positive n, the 
inequality is satisfied when af3 — 1. Therefore, if the 



assertion were always valid, rj would be always zero (or 
very small) for n > 0, which is certainly not the case for 
n > 2 and n < ^ in Fig. [31(e). The assertion is valid 
for the pseudo-Hermitian Dicke model for i < n < 2 of 
which the standard Dicke model corresponding to n = 1 
appears as a special case. 

We conclude with the following: 

(i) Based on our numerical results, we conjecture that 
the level-spacing distribution comes close to the Pois- 
son distribution Pp(s) as the system approaches the 
integrable limit, while for the non-integrable pseudo- 
Hermitian Hamiltonian it should be described by the 
Wigner distribution P\^{s). 

(ii) We have also shown that the assertion ^Td\ that QPT 
is a precursor to a change in the level statistics in the 
standard Dicke model is not valid in general for the 
pseudo-Hermitian Dicke model which has a larger pa- 
rameter space. The assertion holds true for the pseudo- 
Hermitian Dicke model only for a limited range of the 
parameter space. 
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